General theory of measurement with two copies of a quantum state 
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We analyze the possible results of the most general measurement on two copies of a quantum 
state. We show that p can label a set of outcomes of such measurement if and only if there is a 
family of completely co-positive (ccP) maps C M such that the probability of occurrence Prob(p) is 
the fidelity of the map C M , i.e. Prob(p) = Tr(p C M (p)) which must add up to the fully depolarizing 
map. This implies that a POVM on two copies induces a measure on the set of ccP maps (i.e., a 
ccPMVM). We present examples of ccPMVM's and discuss their tomographic applications showing 
that two copies of a state provide an exponential improvement in the efficiency of quantum state 
tomography. This enables the existence of an efficient universal detector. 
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One of the postulates of quantum theory tells us how 
to compute probabilities for the outcomes of measure- 
ments: If the system was prepared in the state p, for 
every outcome p of a measurement there is a projec- 
tor P M such that the probability of occurrence of p is 
the expectation value of the projector in the state 
p, i.e. Prob(p) — Tr{pP ll ). To represent mutually ex- 
clusive outcomes the projectors must be orthogonal and 
they must add up to the identity to ensure that the to- 
tal probability is unity (i.e. P^Py = S^P^, ^ P M 
This postulate, originally formalized by von Neumann 
was extended in the 1970's [2| when the notion of gen- 
eralized measurement was introduced. In such measure- 
ment a positive operator (not necessarily a projector) 
is associated with every outcome p and the probability 
of occurrence of p is Prob(p) = Tr(pA^). The oper- 
ators add up to the identity and define a so-called 
positive operator valued measure (POVM). Neumark's 
Theorem [3j establishes that POVM measurements are 
equivalent to projective measurements for an extended 
system: Any POVM can be implemented via a projec- 
tive measurement on the original system supplemented 
with an appropriately chosen ancillary system. 

In this paper we analyze the predictions of quantum 
theory concerning the results of measurements performed 
when two identically prepared quantum systems are si- 
multaneously available. More precisely, we assume that 



a source produces the state p 



(A,B) 



p ® p ( A and B 



label two systems prepared in the same state p). Our 
goal is twofold: a) To determine the possible distribu- 
tions for the measurement outcomes; b) To present a 
solution to the problem of efficient universal state to- 
mography using copies. We divide the presentation in 
two parts. First we prove the following Theorem (ccP- 
MVM): Given two systems prepared in the same state p 
then p can label a set of possible outcomes of a measure- 
ment on p( A ' B ^ — p p if and only if there is a family 



of completely co-positive (ccP) maps C M such that the 
probability of occurrence Prob(p) is the fidelity of the 
map Cfj,, i.e. Prob(p) — Tr(p (/?)). The maps sat- 
isfy the condition 5Z„ = £ where £ is the map for 
which £(p) = I for any state p. Moreover C p must be 
ccP which means that the composition of C M with the 
transposition must be completely positive. This Theo- 
rem establishes an interesting connection between fami- 
lies of ccP maps and general measurements with copies. 
As a consequence, a measurement with copies defines a 
"ccP map valued measure" (ccPMVM) and vice versa. 
In the second part of this Letter we establish the to- 
mographic power of this type of measurement showing 
that availability of two copies gives an exponential advan- 
tage in solving the general problem of universal quantum 
state tomography enabling us to construct an efficient 
universal quantum state detector and to efficiently esti- 
mate partial purities, concurrences, and other interesting 
quantities. 

To prove the ccPMVM Theorem it is crucial to use the 
so-called Jamiolkowski isomorphism Q that establishes a 
one to one correspondence between linear operators on 
the space H ® H and linear super-operators on the space 
H (a super-operator on H is a map on the space of oper- 
ators over H). For every super-operator C on H we can 
associate an operator C on H (g> H and vice versa. This 
one to one correspondence is realized by the following 
identity: C = (C ® i") (|X) (21) where |2) the unnormal- 
ized maximally entangled state |2) = JV The iso- 
morphism relates positive Hermitian operators on H £g> H 
with completely positive Hermitian super-operators on H . 
In particular, the identity operator is associated to the 
completely depolarizing super-operator £ for which the 
image of every trace one operator is the identity. Using 
this isomorphism the ccPMVM Theorem can be proved 
as follows: An outcome p of a generalized measurement 
on two copies prepared on the state p( A > B ^ = p<g>p is char- 
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Local Bell measurements 
FIG. 1: Proposed scheme for full state tomography. 



acterized by a positive operator C M . The probability of 
such outcome is Prob(p) = Tr(p ® p C^). Jamiolkowski 
isomorphism ensures that for every positive operator 
there is a completely positive super-operator C p such 
that Prob(p) = Tr(p(g>p (C p (g>I)(\l) By replacing 

the explicit form of the state \I) the trace over the sec- 
ond copy can be computed and the probability Prob(p) 
can be rewritten as Prob(p) = Tr(p C ll (p T )), where p T 
denotes the usual transposition of p. Therefore, the prob- 
ability of every outcome of a generalized measurement is 
Prob(p) = Tr{p C M (p)) where C M = C* M o T, with T 
denoting the transposition as a map. Jamiolkowski iso- 
morphism ensures that is completely positive, which 
implies that is completely co-positive. Moreover, as 
the POVM operators add up to the identity, the cor- 
responding super-operators must add up to the com- 
pletely depolarizing map. The relation £oT = T o£ = £ 
completes the prove of the ccPMVM Theorem. 

This Theorem characterizes all measurements with two 
copies and shows that the fidelity of certain families 
of positive super-operators has a direct physical mean- 
ing as it can always be realized as the probability of 
a generalized measurement with copies. We now an- 
alyze an example which is significant from the point 
of view of quantum state tomography. Suppose that 
we have two copies of the state p of an n-qubit sys- 
tem: p( A ' B ) = p ® p. We will assume that we per- 
form a Bell measurement on all pairs formed by the j— 
th qubit of each copy as shown in Figure 1. We write 
the state of each copy as p — J2 qp c q , P T(q,p)/N where 
q = (<7i, . . . , q n ) and p = (pi, . . . ,p n ) are binary n-tuples. 
The generalized Pauli operators T(q,p) are n-fold tensor 
products of the identity and the three Pauli operators 
on each qubit: T(q,p) = X qi Z Pl <g> . . . <g> X q "Z p »(i) qp 
(here qp = y\ qtPk)- Real coefficients c 9iP are such that 
c q>p = Tr(p T(q,p)). 

We can collect the outcomes of all the Bell measure- 
ments in two binary n-tuples (a, b) where (afe,6fe) iden- 
tify the state \f3 ak ,b h ) detected at site k — 1, . . . ,n (Bell 
state \Pa k ,b k ) is an eigenstate of X]~ ® X]~ and Z^ ® Z^ 
with eigenvalue (— l) Qfe and (— l) bk respectively). After 
some manipulation, we can show that the probability 
of occurrence for every possible outcome Prob(a, b) is 
Prob(a,b) = J2q,p(- 1 ) aq+bp+qPc2 q ,p/ N2 ■ These probabil- 
ities are, as the ccPMVM Theorem ensures, the fidelities 



of ccP maps. Indeed, one can show that Prob(a, b) — 
Tr(p C a .b(p)) where the corresponding map is C a ,b{p) = 
T(6, a)p T T(b, a) /N . The ccP character of these maps is 
evident since they are obtained as the composition of the 
transposition with a completely positive super-operator. 
It is useful to analyze the simplest case of a single qubit 
where the coefficients c^o, Co,i and ci ; i are nothing but 
the three Cartesian components of the Bloch vector p 
parametrizing the state as a linear combination of the 
three Pauli's: p = (I + p ■ (?)/2. Then, the maps C a j> 
are such that C a ,b{p) = {I + Pa.b • <?)/4. These opera- 
tors are proportional to states with polarization vectors 
Pa,b = {-l) a p x x + {~l) a+b+1 p y y + (-l)W- Therefore, 
the map C\ t \, corresponding to the singlet |/3i,i), real- 
izes a full inversion on the Bloch sphere. The other Bell 
states have maps corresponding to reflections about the 
three Cartesian planes (where one Cartesian component 
of p changes sign 15(). Adding these four maps we ob- 
tain the fully depolarizing one. Probabilities for the four 
Bell measurements are quadratic in the components of p: 
Prob(a,b) = (l+p -p a ,b)/^ 

Interestingly, the generalized measurement with two 
copies described above can be used to device an efficient 
strategy for quantum state tomography (QST) 0, [|| a 
task whose goal is to extract complete information about 
the state p. This usually requires performing an exponen- 
tially large number of different experiments on identically 



prepared systems. As p = Y^, 



T{q,p)/N is a linear 



combination of A^ 2 generalized Pauli operators T(q,p), 
a straightforward way to determine c 9iP is to measure 
all Pauli's on each qubit. This requires N 2 different ex- 
perimental setups and an exponentially expensive post- 
processing (estimation of each coefficient c q . p typically 
requires the use of a different experimental setup). How- 
ever, we will show that if two copies are available the com- 
plexity of QST is drastically reduced. The use of copies 
was analyzed before in this context and it was shown 
that it enables us to extract more information than single 
copy measurements @, BSLIS U, 13 • Here we general- 
ize these results showing that the number of experiments 
required to gather the information to estimate every coef- 
ficient c 2 with fixed precision is independent of the num- 
ber of qubits. A possible strategy to determine all c q p 
is to estimate all probabilities Prob(a, b) and later to ob- 
tain such coefficients by performing a Hadamard trans- 
form. However, this would be highly inefficient method 
since to obtain such coefficients we would need to es- 
timate an exponentially large number of exponentially 
small probabilities (it is simple to see that probabilities 
obey Prob(a,b) < l/N). On the contrary, to efficiently 
estimate c 2 , p we can proceed as follows: After perform- 
ing the Bell measurements described in Figure I we can 
multiply the detected values of the operators Xj <g> X, 
and Zj (g) Zj to obtain the value of any Pauli operator of 
the form T(q,p) ® T(q,p). By repeating this experiment 
a number of times we can estimate the corresponding ex- 
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pectation value of such operator, which is precisely what 
we need to compute \c q . p \ since: 



(T(q,p)®T(q,p)) p „ =c 



<?,p' 



(1) 



We can show that the number of experimental runs, Me, 
required to estimate any c q:P (up to a sign) with fixed 
precision is independent of the number of qubits and 
is only fixed by the precision. In fact, every measure- 
ment yields binary values for T(q,p) <8> T(q,p). There- 
fore, after Me repetitions we compute the average re- 

The central limit theorem im- 



sult that we denote c qp 



for 



satisfies 



plies that the standard deviation a qiP . ,, 

&q,p < l/V Me- Thus, there is a number k such that 
Cg p € [c 2 p — ka q ^ p ; c q p + ka q ^ p ] with probability p. This 
bound propagates to |c 9jP | that with the same probabil- 
ity p will be found in an interval centered at |c gjP | with 

a width i~ q ' P i ■ On the other hand, if one wants to 

2 \ c q,p\ 

estimate each \c qtP \ larger than a fixed S with an uncer- 
tainty e, and obtain a correct value with probability p, 
the number of required repetitions is: Me > k 2 /A5 2 e 2 
where k is chosen to satisfy p = ei ^{^^j- Thus, the 
number of repetitions does not depend on n but only on 
the precision e, the minumum measurable value 5 and the 
probability of success p. The method is "quantum effi- 
cient" as the number of quantum resources (i.e., copies of 
the quantum state, measurements, etc) is constant given 
a required precision. However, classical resources to de- 
termine every c q , p are still exponential in n due to the 
fact that there are N 2 such coefficients. 

The setup of Figure 1 is a universal quantum state de- 
tector that, as opposed to previously proposed ones, is 
efficient. Using it to estimate a set of coefficient with 
fixed precision requires only a number of experiments 
(and classical resources) which is independent of n. Uni- 



versal state detectors were introduced before [13j but, as 
we will show now, they are inefficient. These universal 
detectors do not use copies but ancillary systems pre- 
pared in a known state p = ^}2 q p c q °) p T(q 1 p) / N . When 
systems (A) and (B) are respectively prepared in states 
p and po, it is simple to show that by performing joint 
Bell measurements on every pair of qubits we obtain 



- c (0) 



(T(q,p)®T(q,p)) 



Therefore, knowing c q °) p and measuring the expectation 
value appearing in ^ we can determine c q p . How- 
ever, a universal detector must use a state po that is 
not related with p. This is the origin of the ineffi- 
ciency of the method, as can be seen as follows. Clearly, 
we can determine c q/p only if the corresponding c q ° p is 
non- vanishing. Moreover, the smaller the value of c q °) p 
the higher the precision required in the estimation of 



(2) 



,(o) 



are maximal. This is the case for stabilizer states 



(common eigenstates of a commuting set of TV Pauli oper- 
ators T(q,p)s). For such state there are N non- vanishing 
coefficients c q °l s taking values equal to ±1. For such p 
the universal detector can only be used to estimate N 
c 9>p' s providing no information about the N 2 — N remain- 
ing ones, denying its universality. On the other hand all 
the coefficients could be estimated using a state po with 
non- vanishing c q °l for all (q,p)'s. The problem for such 
unbiased po is that all <4°p's are exponentially small. The 
reason for this is that T, q , p c(0 \,p/ N < Therefore, 

each coefficient c q °) p is 0(1/ y/N). Then, if we use © 
to estimate them with fixed precision we need exponen- 
tially high precision in the estimation of the expectation 
value (T(q,p) ® T(q,p)). For this reason the method is 
inefficient (the universal detector would have to be used 
an exponentially large number of times to achieve a fixed 
precision) . Clearly, the use of a copy instead of an ancilla 
provides a simple way out of this problem. 

Full quantum state tomography is always exponen- 
tially hard as the number of unknown parameters scales 
as 4". Therefore, it is crucial to conceive efficient meth- 
ods for partial characterization of quantum states. Re- 
markably, the strategy described in Figure I is effi- 
cient also for this purpose. To see this we consider 
"coarse grainded" Bell measurements: For any Bell 
state \Pm,n) we can estimate the probability to de- 
tect an even (odd) number of them in the measure- 
ment of all pairs. Then, we can compute AProb m ^ n = 
Prob(even # |/3 m , n )) - Pro6(odd # |/3 m ,„)). It is simple 
to show that AProb mn — Tr(O m>n (p) p) where the (not 
necessarily positive) map O m , n is such that 



AProb m . n = 



N ^ 



(3) 



Here, the N 2 components of the vector s q ^ p are s q . p = 

(_ 1 )(m+l)(a a ,+a I( )(_ 1 j(n+l)(a >[ +a y ) j where ^ ( reS p.' a y , 

a z ) denotes the number of qubits for which the Pauli 
operator T(q,p) contains an X (resp. Y, Z) operator. 
For example, for the singlet |/3n), s 9iP = 1 for every a. 
For any other Bell state half of the components of s 9jP 
are equal to +1 and the other half are equal to —1. For 
the singlet the above formula reduces to 



APr 6 M = ^E c L = Tl '(^ 

9:P 



(4) 



(T(q,p)®T(q,p)) 



p®Po ' 



Consider first a state for which 



Thus, this measurement reveals the purity of the state. 
Partial purities can be detected in the same way: Con- 
sider the state pj, obtained after tracing out the qubits 
for which the binary n-tuple J is zero. Purity of such 
state is the sum of c 2 p for the coefficients associated with 
Pauli operators containing the identity in the qubits for 
which the corresponding component of J is zero. To ob- 
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tain it we must use ((4]) counting singlets only in the qubits 
where the corresponding bit of J is equal to unity. 

The above method for estimating purity is equiva- 
lent to the one proposed by Ekert et al who used the 
fact that purity is equal to the expectation value of the 
swap operator in the state p^ A ' ® p^ [14| . But our 
results also show that by making more general coarse 
grained Bell measurements we are not only able to ef- 
ficiently detect partial purities. In fact, with the same 
effort we reveal other quantities that partially charac- 
terize the quantum state and have the form AP5 = 

(E {q , p) es c lp-^{ q ,p)€S c ip) l N > where i S >S} is a 
partition of the iV 2 coefficients c 9iP in two halves. It is 
possible to generalize this even further by grouping Bell 
states in each pair of qubits (in this case one can attain 
linear combinations with vectors s q „ that have a differ- 
ent number of ±1 components). Other wheighted sums 
of squares of c a over certain sets of (q,p)'s can also be 
obtained in this way. 

It is interesting to consider another related coarse 
grained Bell measurement: we can estimate the prob- 
ability Pm,n to find all pairs of qubits in the subspace 
orthogonal to |/3 mj „). For the case of the singlet this is 
related to the multipartite concurrence that for a pure 
n-qubit state is [ll| 



C(p)=2 1 -/y(2"-2)-^Tr(p 2 ). (5) 

Here the n-tuple I labels every nontrivial subset of the n 
qubits. As noticed in this can be rewritten as: 

C(p) = 2yfl -p[ a } l \ (6) 

therefore, concurrence of a pure state can be esti- 
mated efficiently using the type of coarse grained Bell 
measurements we described above. The probability 
Pii^ is also a quadratic form of the coefficients c 9iP : 

_ Y^, qp c 2 p i a °/N 2 , where ao is the number of 
qubits for which T(q,p) contains a factor equal to the 
identity. More generally, if we measure Pm U n we ob- 
tain a quantity that is not a concurrence but pro- 
vides different tomographic information. It can be ex- 
pressed as p)n$ = ^q^c^pfq^p/N 2 ! where f 9iP = 

As we mentioned, for an efficient partial QST it is 
necessary to estimate certain differences between prob- 
abilities such as in (3). Such probabilities, as the ccP- 
MVM Theorem states, are fidelities of ccP maps. Thus, 
the right hand side of equation (3) can be expressed as 
the fidelity of a map which is the difference between two 
ccP maps: AP m>n = Tr(pO m , n (p)) (O m , n is not ccP). 
One of such maps is the identity, which is not ccP since 
when composed with the transposition becomes non-cP. 
For that case, one finds AP1.1 = Tr(p 2 ). Thus, purity 



can never be evaluated as the probability of a generalized 
measurement but only as the difference between probabil- 
ities. In turn, fidelities of positive maps that are not ccP 
can be related via the Jamiolkowski isomorphism with 
expectation values of Hcrmitian operators (not necessar- 
ily positive). Such Hcrmitian operator can be written as 
the difference between two positive operators. For the 
example disc ussed above, where the positive map is the 
identity, it turns out that the non-positive operator is the 
one implementing the swap between the two copies as the 
well known identity AP^i = Tr{p 2 ) = Tr(Swap p ® p) 
shows. Summarizing, we showed that by performing gen- 
eral measurements on two copies of a quantum state one 
always detects probabilities that are fidelities of com- 
pletely co-positive maps. We also described techniques 
to perform efficient quantum state tomography by taking 
advantage of measurements with copies. 
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